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Abstract 

We have used the Planck all-sky submillimetre and millimetre maps to search for rare sources distinguished by extreme brightness, a few hundred 
millijanskies, and their potential for being situated at high redshift. These “cold” Planck sources, selected using the High Frequency Instrument 
(HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame 
far-infrared peak redshifted to the frequency range 353 - 857 GHz. This colour-selection favours galaxies in the redshift range z = 2^1, which 
we consider as cold peaks in the cosmic infrared background. With a 415 beam at the four highest frequencies, our sample is expected to include 
overdensities of galaxies in groups or clusters, lensed galaxies, and chance line-of-sight projections. We perform a dedicated Herschel- SPIRE 
follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500/im sources, in comparison to reference SPIRE fields. About 
94 % of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350 pm, with 3 % peaking at 500 pm, 
and none peaking at 250 pm. About 3 % are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at 
redshifts z > 2.2. Only four targets are Galactic cirrus, yielding a success rate in our search strategy for identifying extragalactic sources within the 
Planck beam of better than 98 %. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance 
above 10 a. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z - 2, assuming a single common dust temperature 
for the sources of Ta = 35 K. Under this assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4 x 10 12 L 0 , yielding 
star formation rates of typically 700 M 0 yr~‘. If the observed overdensities are actual gravitationally-bound structures, the total IR luminosity of 
all their SPIRE-detected sources peaks at 4 x 10 13 L Q , leading to total star formation rates of perhaps 7 x 10 3 M 0 yr~' per overdensity. Taken 
together, these sources show the signatures of high-z (z > 2) protoclusters of intensively star-forming galaxies. All these observations confirm the 
uniqueness of our sample compared to reference samples and demonstrate the ability of the all-sky Planck- HFI cold sources to select populations 
of cosmological and astrophysical interest for structure formation studies. 

Key words. Galaxies: high-redshift, clusters, evolution, star formation - Cosmology: observations, large-scale structure of Universe - 
Submillimetre: galaxies - Gravitational lensing: strong 


1. Introduction 

_ Cosmological structure formation in the linear regime is now 

* corresponding author e-mail: herve. dole@ias. u-psud. fr fairly well constrained by observations at the largest scales (e.g., 
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iTegmark et al.ll2004l ICole et al .1120051) . In the non-linear regime 
the situation is less clear because structure formation becomes a 
complex interplay between dark matter collapse and the hydro¬ 
dynamics of baryonic cooling. In the particularly dense environs 
of the most massive dark matter halos, this interplay should lead 
to vigorous episodes of rapid star formation and galaxy growth, 
giving rise to copious amounts of FIR (far-infrared) emission 
from dust heated by young stellar populations in massive galax¬ 
ies during periods of intense star formation. 

From the standpoint of galaxy evolution, studying this 
intense star formation epoch in massive dark matter halos 
may provide a wealth of observational constraints on the 
kinematics and evolutionary history of galaxies in massive 
galaxy clusters, as well as the intracluster gas. From the 
point of view of cosmology, clusters yield information on 
non-Gaussianities of primordial fluctuations and can challenge 


the ACDM (cold dark matter) mo d el dBro dwin et al. 
Hutsil 20lol Williamson et al.l 1201 it Harrison & Colesl 


Holz & Perlmutter 2012t Waizmann et al.112012 : iTrindade et al.l 

20131) . while lensed sources act as probes of dark matter ha¬ 


los, and both are probes of cosmological parameters such as 
Qm (matter density today divided by the critical density of the 
Universe) and cr 8 (the rms fluctuation in total matter - baryons 
+ CDM + massive neutrin os - i n 8 h 'Mpc spheres today at 
z = 0) JPlanck Collaboration XVll 120141) . Furthermore, clus¬ 
ters of galaxies are crucial objects that bridge astrophysics and 


ing the measurement of <x« (Planck Collaboration XVII 

2014 

Planck Collaboration XX 

2014 

Planck Collaboration I 

2015 

Planck Collaboration XIII 

2015 

), all of which has led to a de- 


bate between cluster phenomenology and cosmological physics. 

The extragalactic sky in the submillimetre (submm) and 
millimetre (mm) regime has been of considerable scientific 
interest for over two decades, with the distinct advantage 
that the steep rise in the redshifted Rayleigh-Ieans tail of 
the modified blackbody emitted by the warm dust in in¬ 
frared galaxies largely compensate s for cosmologica l dim¬ 
ming, the “negative k-correction” ( France schini et afl 1199 It 


iBlain & Longairl 119931 : iGuiderdoni et al.i 11997h . As a conse¬ 
quence, the flux density of galaxies depends only weakly on 
redshift, opening up a particularly interesting window into 
the high-redshift Universe (typically 2 < z < 6). Constant 
improvements in bolometer technology have led to impres- 


based, balloon and space-borne telescopes (e.g., Hughes et al. 

1998; 

Barger et al.| 11998; Chapman et al.l 

2005: lLagache et al. 

2005; 

Patanchon et al. 2009; Devlin et al 

2009; Chapin et al. 

2009; 

Negrello et al. 2010; Vieira et al l 2010; Oliver et al.l 2010; 

Eales et all2010|). 


Nevertheless, only with the recent advent of wide-field sur¬ 
veys with astrophysical and cosmological scope have the sys¬ 
tematic searches become efficient enough to identify the bright 
est of these objects with flux densities above about lOOmJy 
at 350 pm, e.g., with He rscheH\ dPilbratt et alj|2010t lEales et al 


201L; lOliver et ~aT1l2Q 1 Ot ), the South Pole Telescope ( Vieira et al 
20iol) , WISE (I Wright et all f20Tok [Stanford et al.l 12014 . and 


Spitzer dPapovichl 120081 : IStanford et alJ l2012h . Such sources 
are very rare. For example, the surface density of red sources 
brighter than 300 mJy at 500 pm is 10 2 deg 2 for strongly 
lensed galaxies, 3 x 10 2 deg 2 for AGN (active galactic nu¬ 


1 Herschel is an ESA space observatory with science instruments 

provided by European-led Principal Investigator consortia and with im¬ 

portant participation from NASA. 


clei, here radio-loud, mostly blazars), and 10 1 * * d eg 2 for late- 


type galaxies at moderate redshifts (e.g., I Negrello et al 


20101). Other models predict similar tren ds ( Paciga et al 
Lima et a~fll2010t iBethermin et al.l 1201 it iHezaveh et alJ 


200 ' 


200 ' 

20121 ) . 


This makes even relatively shallow submm surveys interesting 
for searches of high-redshift objects, as long as they cover large 
parts of the sky. Studies of gravitationally lensed g alaxies a t high 
redshifts originating from these surveys dNegrello et all 201C; 

Combgs^alj |2()l4^ Bussmann et al ]|2013nHerranz et aD 2013 


Rawle et al1l2014t [Ward low et al.ll2013l) illustrate the scientific 


potential of such surveys for identifying particularly interesting 
targets for a subsequent detailed characterization of high-z star 
formation through multi-wavelength follow-up observations. 

The power of wide-field surveys in detecting the rarest ob¬ 
jects on the submm sky is even surpassed with genuine all¬ 
sky surveys, which systematically and exhaustively probe the 
brightest objects in their wavelength domain down to their com¬ 
pleteness limits. Here we present a search for the rarest, most 
extreme high-redshift (z, >2) candidates on the submillime¬ 
tre sky , which was p erformed with the Plancl Q all-sky survey 
dPlanck Collaboration lll201~4l 20151). Th e Pl anck Cata logue of 


Compact Sources (PCCS. lPlanck Collaboration XXV11 fl l2014l) 

has a completeness limit of about 600 mJy at the highest fre¬ 
quencies, which corresponds to Lfi r - 5 X 10 13 L Q at z — 2. 
With a 5' beam (iPlanck Collaboration VIll l2014t) at the four 
highest frequencies (corresponding to a physical distance of 
about 2.5 Mpc at z — 2), we expect that sources with bona 
fide colours of high -z galaxies in the Planck -HFI bands are ei¬ 
ther strongly gravitationally lensed galaxies, or the combined 
dust emission of multiple galaxies in a shared vigorously star¬ 
forming environment in the high redshift Universe. The lat¬ 
ter case is consistent with the result that submm galaxies or 


tered 

(Blain et al. 2004; Farrah et al. 2006; 

Magliocchetti et al. 

2007; 

Austermann et al. 2009; Santos et al. 

2011; Ivison et al. 

2012; 

Valtchanov et al. 2013; Noble et al. 201J; Clements et al. 


20141) . and may include massive galaxy clusters during their ma¬ 


jor growth phase. It is also possible of course that the sources are 
chance alignments of multip le high-redshift galaxies projected 
onto_the same line of sight dNegrello et al.ll2005l 120071120101 : 
Chiang et al. 2013), or of multiple, lower-mass galaxy groups or 
clusters. 


Identifying high-redshift cluster candidates directly by the 
signatures of their total star formation is a very useful com¬ 
plement to the diagnostics used to identify galaxy clusters 
so far. Most systematic searches today focus on the pri¬ 
mary constituents of more evolved, lower-redshift clusters, 
like their populations of massive, passively evolving galax¬ 
ies (’red-sequence galaxies’), the hot intracluster medium (ei¬ 
ther through X-ray emission or the Sunyaev-Zeldovich ef¬ 
fect), or the suspected progenitors of today’s bright est c lus- 


ter galaxies, in particular high-redshift galaxies (Steidel et al. 

2000 

Brand et al. 2003 

Kodama et al. 2007; Scoville et al. 

2007 

Venemans et al.l 

2007; Papovich 2008t Daddi et al. 

2009 

Brodwin et al.l 2010; Galametz et alj2010; Papovich et al. 

2010 

Capak et al. 

201 It Hatch etal. 2012 Kuiperetal. 2011; 

Ivison et al.l l2013t 

Muzzin et al. 2013; Wvlezalek et al. 2013a; 


2 Planck is a project of the European Space Agency - ESA - with in¬ 
struments provided by two scientific Consortia funded by ESA member 
states (in particular the lead countries: France and Italy) with contri¬ 
butions from NASA (USA), and telescope reflectors provided in a col¬ 
laboration between ESA and a scientific Consortium led and funded by 
Denmark. 
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Chiang et al.1 20141: Cooke et al 


2014 ICucciati et all [2014;, 


Mei et all 2014 iRettura et al.1 2014 . for instance). With the lim¬ 
ited sensitivity and the large beam of Planck, in turn, we effec¬ 
tively select the most intensely star-forming Mpc-scale environ¬ 
ments in the high-redshift Universe. 


Using Planck, we have selected putative high-redshift ob¬ 
jects with spectral energy distributions (SEDs) of warm dust, 
peaking between observed frequencies of 353 and 857 GHz. 
This has a net effect of selecting sources either peaking at 
545 GHz, or to a lesser extent sources having their infrared 
peak between 353 and 857 GHz. This equates, in principle, 
to redshifted infrared galaxies at z ~ 2-4. We refer to these 
as “cold” sources of the Cosmic Infrared Background (or 
CIB), wh ich have a red and thus potentially redshifted SEP, 
The CIB dPuget et ahlfl996l:lHauser et ahlfl99llHauser & Dwekl 


|200T1 [PoU et al.ll2004]pianck Collaboration XXXll2014l) is the 
integrated relic emission in the broad infrared range, typi¬ 
cally 8 gm to 1 m m, where the emission reaches a maximum 
(Dole et al.1120061) . Physically, such objects correspond to galaxy 
and AGN formation and evolutionary processes, and more gen¬ 
erally the history of energy production in the post- rec omb i natio n 
Universe (e.g., iKashlinskvl 120051 : iPlanck Collaboration XVIlH 
1201 lh . The CIB, as observed in the submillimetre, is consid¬ 
ered _as_a_proxy for intense star formation at redshifts z > 1 
(Planck Collaborat ion XXXll 2014 as well as for the mass of 
structures (Planck Collaboration XVlifll20l4 . 


The Planck c oll aborati on has also released the PCCS 
(Planck Collaboration XXVHll20l4 . containing 24381 sources 
at 857 GHz, with 7531 sources at Galactic latitudes \b\ > 30°, of 
which many are of interest for extragalactic studies. The main 
difference between the cold sources of the CIB and the red 
sources of the PCCS (both used in this work) can be summarized 
as a threshold difference: cold CIB sources are detected in the 
CIB fluc tuati ons with a com bined spectral and angular filtering 
method ( Montier et al.ll20f(X ). while PCCS sources are detected 
independently in the frequency maps u sing an angular filtering 
method with a higher threshold (Planck Collaboration XXVHll 
12014 . 


This paper presents the observationsjind anal ysis of our ex¬ 
tensive dedicated Herschel -SPIRE (Griffin et al.11201(f) follow¬ 
up of 234 Planck sources (either selected from the CIB fluctu¬ 
ations or from the PCCS). The paper is structured into seven 
sections. In Sect. [2] we detail the Planck parent sample and 
the Herschel observations. Section [3] gives a technical descrip¬ 
tion of the algorithms used in the generation of the SPIRE 
photometry and the catalogue. In Sect. [4] we use statistics to 
characterize the Herschel observations; in particular we quan¬ 
tify the overdensities and the colours of the SPIRE counterparts 
and propose a classification of either overdensities or lensed 
candidates. In Sect. 0 we discuss the properties of the lensed 
source candidates, while in Sect. [ 6 ] we focus on the overden¬ 
sities and their characterization, including a stacking analy¬ 
sis. Conclusions are reported in Sect. |7] The Appendices con¬ 
tain information on the SPIRE catalogue generation and the 
number counts, and a gallery o f sample fields. We use the 
Planck 2013 cosmology (Planck Collaboration XVIl 12014 ta¬ 
ble 5: P/anck+WP+highL+BAO) throughout the paper. 


2. Sample selection and observations 

2.1. Planck observations and selection 

Planck 1^] observ ed the whol e sky at frequencies between 30 and 
857 GHz (Planck Collaboration 12014 ). We made two different 
selections to follow-up with Herschel. first, using the maps and 
looking for cold sources of the CIB; second, using the public 
catalogue of compact sources (PCCS). 


2.1.1. Cold sources of the CIB in the Planck maps 


We mak e use o f the Planck -HFI (High Frequency Instrument, 
[Planck Collaboration VI 2014 data as well as the /PAS'/IRIS 
data ( Miville-Deschenes & Lagachell2005l) at 100/mi. For this 
purpose we use the cleanest 26% of the sky (in the Planck 
857 GHz map), defined by a minimal cirrus contamination, 
/V H i < 3 x 10 2( > cm 2 . The detection and selection algorithm, 
based on iMontiere t al.l (20 1 Oh can he summa rized in the fol¬ 
lowing seven steps ( Planck Collaborationll2015l) . 


(i) CMB cleaning: the 143 GHz Planck-HPl map is extrapo¬ 
lated to the other bands according to a CMB spectrum and 
removed from the maps at other frequencies. 

(ii) Galactic cirrus cleaning: the IRAS/IRIS 100/rm map, con¬ 
sidered as a “warm template” of Galactic cirrus, is extrapo¬ 
lated to the Planck -HFI bands, taking into account the local 
colour, and removed from the maps following the prescrip¬ 
tion of the CoCoCoDeT algorithm (Montier et al.llioiol) . 

(iii) Construction of excess maps: the 545 GHz excess map is 
defined as the difference between the cleaned 545 GHz 
map and a power law interpolated between the cleaned 
maps at 857 GHz and 353 GHz. 

(iv) Point source detection in the excess maps: we 
apply a Mexican hat_ type detection algorithm 
(Gonzalez-Nuevo et al.l 120061) with a size parameter 
of R — 5' in the 545 GHz excess maps. 

(v) Single frequency detection: we also require a detection in 
each cleaned map at 857 and 353 GHz. 

(vi) Colour-colour selection: we apply two criteria on the 
‘S ’545 /•S '857 (he., 545 GHz to 857 GHz) and S 353/£545 (i.e., 
353 GHz to 545 GHz) flux density ratios to select the red¬ 
dest sources. 


This produces a dozen hundred candidates on the cleanest 
26 % region of the submm sky. We note that this procedure has 
been set-up early on in the Planck project. The final catalog 
of Planck high-z candidates is being generated using a similar 
(but not exactly iden tical) method which will b e described in 
a forthcoming paper (Planck Collaborationll2015l) . In particular, 
the CMB estimate will not be the 143 GHz HFI map, but instead 
the CMB derived by component separation. The present paper 
focuses on the first candidates followed up by Herschel. 


2.1.2. Planck PCCS sources 

We make use of the PCCS to choose a sample of high-z 
sources selected by the expected peak in their thermal dust 
spectrum in the rest-frame far infrared. Our four step proce¬ 
dure here is based on the work of iNegrello et al j (2010i) . First 
of all, we use the 857 GHz Galactic mask, keeping 52 % of 

3 Planck data (maps and catalogues) can be 
downloaded from the Planck Legacy Archive 
http://pla.esac.esa.int/pla/aio/planckProducts.html 
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Figure 1. SPIRE maps at 250, 350, and 500 pm of one typical overdensity in our sample (a Planck cold source of the CIB - Cosmic 
Infrared Background). The thick white contour at 350 pm shows the Planck IN region at 857 GHz, and the same is true for the map 
at 500 pm at 545 GHz (i.e., the 50 % contour of the Planck source maximum). Thin contours correspond to the overdensity of SPIRE 
sources at each wavelength, marked at 2cr, 3cr, etc. (see Sect. I4.2l for details). SPIRE identifies a few (typically 5 to 10) sources 
inside the Planck contour. We use these contours to separate the IN and OUT zones. These data come from a 7' x 7' SPIRE scan 
(see Tableland Sect. 12.3k 


PICK 02448*549 


a 


a 


,1 


a 


Figure 2. SPIRE maps at 250, 350, and 500yum for a lensed source. In this case a single very bright SPIRE source is detected. This 
source is co nfirmed at the Jy level and has a redshift of z = 3.0, measured from multi-line spectroscopy acquired at IRAM using 
EMIR (ICanamerasll2015l) . The white contour is the Planck IN region, defined in the same way as in Fig. Q] These data come from a 
7' x 7' SPIRE scan (see TableQ]and Sect. 12.3k 



350,4m | 


250^m 


500,«m 


the sky dPlanck Collaboration Int. Vlll2013h . Secondly, we se¬ 
lect all the sources with S/N > 4 at 545 GHz and the colours 
S 857 /.S’ 545 < 1.5 (where S v is the flux density at frequency 
v in GHz) and 5 217 < 5 353 . Thirdly, we inspect each source 
with respect to the NASA/IPAC Extra galactic Database (NED), 
IRAS maps dNeugebauer et al.l i 19841) . and optical maps using 
ALADIN. Any object identified as a local galaxy, a bright radio 
source or Galactic cirrus is removed (about half of the objects). 
Finally, we remove any PCCS source already falling in the H- 
ATLAS or South Pole Telescope (SPT) survey fields. 


2.2. The Planck parent sample and three subsamples 

Planck and Herschel were launched on 14 May 2009, and 
started routine scientific observations a few months later 
(tPlanck Collaboration 1)12014 1 IPilbratt et al.H2010 ). Our Planck 
parent sample has evolved with time as a function of the de¬ 
livery of internal releases of the intensity maps (typically ev¬ 
ery 6 months). We thus built three samples: one for each 
main Herschel Open Time call for proposals (OT1 and OT2), 
and a last one for the Herschel ’must-do’ DDT (Director’s 
Discretionary Time) observations (which we call HPASSS for 
Herschel and Planck All Sky Legacy Survey Snapshot). We re¬ 
quested only SPIRE data, as it is better optimized than PACS for 
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Table 1. SPIRE programmes following-up Planck cold sources 
of the CIB. HPASSS is composed of 124 sources selected from 
the maps (two were repeated, so 126 fields), together with 28 
from the PCCS (four were cirrus dominated, so 24 net fields). A 
total of 228 sources are useable (Sect. lT2t . 


Programme No. of fields Map size AOR 11 time 


OT1 . 10 15' x 15' 3147 s 

OT2 . 70 7' x 7' 838 s 

HPASSS... 124 (126) (9 b ) + 24 (28) 7'x 7' 838 s 

Total. 228 (234) 


a Astronomical Observation Request. 
b Nine HPASSS sources are from archival data. 


the targeted redshift range. Our main goal was to identify the 
counterparts contributing to the Planck- HFI detections, hence 
the focus on rapid SPIRE follow-up of a maximum number of 
targets. Obtaining PACS observations would have enhanced the 
angular resolution and wavelength coverage, but for these short 
observations would not have resulted in further constraints in 
most cases. 

The samples are summarized in Table [Q together with the 
map sizes and time per Astronomical Observation Request 
(AOR). The OT1 sample refers to the Herschel first call 
for Open Time observations in July 2010 (P.I.: L. Montier). 
The OT2 sample refers to the Herschel second call for 
Open Time observations in September 2011 (PI.: H. Dole). 
Figs. [I] and [ 2 ] show examples of this observation strategy. 
The Herschel and Planck All-Sky Source Snapshot Legacy 
Survey (HPASSS) was set up by the Planck collaboration (P.I.: 
H. Dole) in the response to ’must-do’ DDT programme in 
June 2012. At that time, the preliminary PCCS was internally 
released to the collaboration, and we could benefit from the 
unique and timely Herschel and Planck synergy when the 
Planck products were approaching their final state and when 
Herschel surveys already demonstrated the efficiency of the 
SPIRE observations. We note that the SPIRE data from 9 
Planck fields included in the HPASSS sample are Herschel 
archival data: 4C24.28-1 (program: OT2_rhuub22); NGP_vl, 
NGPJil, NGP_h2, SGP_sm3-h, GAMA12_rnl, NGP_v 8 , 
NGPJ 16 (KPOT_seales01 _2); Lockman_SWIRE.offset.l-l 
(SDP_soliver_3); and Spider-1 - (OT1 _mmiville_2). 

2.3. SPIRE data processing, total SPIRE sample and 
definition of IN and OUT Planck regions 

The SPIRE data were processed starting with ’Level 0’ using 
HIPE 10.0 (Oil; 2010 ) and the calibration tree ’spire.cal.lO.l’ 
using the Destriper modukQhttps://nhscsci. ipac.caltech.edu as 
mapmaker (baselines are removed thanks to an optimum fit be¬ 
tween all timelines) for most observations. For 33 AORs, we 
had also to remove some particularly noisy detectors (PSWB5, 
PSWF 8 , and PSWE9 affected for Operational Days 1304 and 
1305) from the Level 1 timelines. In that case, we processed the 
data with the naive scan mapper using a median baseline removal 
(the destriper module worked with all bolometers in Level 1). 
Turnarounds have been taken into account in the processed data. 
The useable sky surface area is thus extended, and goes beyond 
the nominal 7'x7' or 15'x 15' as specified in the AORs (Table|T}. 


Table 2. SPIRE 1 cr total noise (instrument + confusion) levels 
measured in various fields, in mJy. Our Planck fields are denoted 
’Planck high-z’. 


1 cr noise level 


Field 

250 pm 

350 pm 

500 pm 

Planck high-z. 

9.9 

9.3 

10.7 

HerMES Lockman-SWIRE_ 

10.1 

10.5 

12.0 

HLS. 

14.1 

12.6 

14.2 


Since we are close to the confusion limit, we included a check 
to confirm that the non-uniform coverage imposed by including 
the edges does not change the source detection statistics. 

We have a total of 234 SPIRE targets. Of these, two fields 
were repeated observations, and have been used to check the ro¬ 
bustness of our detections. Four fields from the PCCS appear as 
cirrus structures: diffuse, extended submillimetre emission with¬ 
out noticeable point sources. These fields were removed from 
the sample. This means that four out of 234 fields (1.7 %) were 
contaminated by Galactic cirrus. The success rate of avoiding 
Galactic cirrus features is thus larger than 98 %, thanks to our 
careful selection on the Planck maps of the cleanest 35 % of the 
sky. 

Our final sample thus contains 228 fields (i.e., 234 minus 
two repeated fields minus four cirrus-dominated fields). They 
are composed of (See TableQ}: 10 sources from OT1; 70 sources 
from OT2; 124 objects from HPASSS CIB; and 24 sources from 
HPASSS PCCS. 

Each SPIRE field of a Planck target is then separated into 
two zones: IN and OUT of the Planck source at 545 GHz (the 
frequency where our selection brings the best S/N ratio). The IN 
region boundary is defined as the 50 % Planck intensity contour, 
i.e., the isocontour corresponding to 50% of the peak Planck 
flux, and encompasses the Planck-HFl beam. The OUT region 
is defined to be outside this region (see the thick white contours 
in Figs.0and[2]l. 

2.4. Ancillary SPIRE data sets 

For a first characterization of our sample, and to infer whether it 
is different from other samples of distant galaxies observed with 
SPIRE, we will compare with the data obtained in other SPIRE 
programmes. Since we suspect that most of our targets contain 
overdensities like proto-clusters of galaxies, it will be useful to 
contrast it with samples of galaxy clusters at lower redshift, as 
well as against blank fields. Our two comparison samples are as 
follows. 

1. pie HerMESB (lOliver et alJ 1201 Oh and H-ATLAS0 
(lEales et al.l l2010h public data as reference fields. In 
particular, we will be using the ’level 5’ Lockman HerMES 
field, which has a similar depth to our SPIRE observations. 

2. The SPIRE s napshot programm e of local or massive 
galaxy clusters (lEgami et al.ll2010l) including: the “Herschel 
Lensing Survey” (HLS); the “SPIRE Snapshot Survey of 
Massive Galaxy Clusters” (OT1); and the “SPIRE Snapshot 


5 http://hedam.oamp.fr/HerMES/ 
Scan map destriper details: https://nhscsci.ipac.caltech.edu/sc/index.php/Sphtffhpt^BHnrMapiidsliBpoirg/ 
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250 /<m 350/un 500/un 



100 1000 100 1000 100 1000 

^input 250 (mJy) 

^input 350 (mJy) 

^input 500 (mJy) 


Figure 3. Photometric accuracy from Monte Carlo source injections in real maps at 250, 350 and 500gm. The 3 cr (dotted line) and 
4 cr (dashed line) levels are also shown. See Sect. l3.2l for details. 


500/un 350 250 /im 



S 500 [mJy] s 350 [mJy] S 250 [mJy] 


Figure 4. Completeness levels from Monte Carlo source injections in SPIRE maps of different data sets at 250, 350 and 500 /am. 
Plotted are our Planck fields (black), HerMES sources (orange), and Herschel Lensing Survey (HLS) clusters (blue). The 50% 
(dotted line) and 80% (dashed line) completeness levels are also shown. See Sect. l3.2l for details. 


Survey II, Using SPT/CODEX Massive Clusters as Powerful 
Gravitational Lenses” (OT2fl 

Finally, we use the OT2 data from “The highest redshift 
strongly lensed dusty star-forming galaxies” (PI. J. Vieira) pro¬ 
gramme to carry out technical checks, since those AORs are very 
similar to ours. 

3. SPIRE Photometry 

3.1. Photometric analysis 

We first measure the noise in each map at each wavelength 
(channel) by fitting a Gaussian function to the histogram of the 
maps and derivin g the standard deviation, gyhannd, a m ixture of 
confusion noise dDole et al.ll200.3i : iNguven et al.ll2010h and de¬ 
tector noise. The values are reported in Table[2] Data are sky sub¬ 
tracted and at mean 0, as done in HIPE 10. All pixels are used to 
plot the noise histogram. Histogram have a Gaussian shape, with 
a bright pixel tail. The bright pixel tail contribution (larger than 3 

7 Herschel public data can be down¬ 
loaded from the Herschel Science Archive: 
http://herschel.esac.esa.int/Science_Archive.shtml 


sigma) to the histogram is nearly 1 % in number for each field at 
350um. Then we extract the sources using a blind method and by 
a band-merging procedure described in the following sections. 

3.2. Blind catalogues 

We detect blindly, i.e., independently at each frequency, 
the sources using StarFinder dDiolaiti et al.l 120001) . We use 
Gaussian point spread functions (PSFs) with FWHM of 18.15" 
for 250/rm, 25.15" for 350yum, and 36.3" for 500 /m0. These 
individual band catalogues are used for checking photometric 
accuracy and completeness, and to produce number count esti¬ 
mates. The catalogues are then band-merged in order to quantify 
the colours of the sources. 

We employ Monte Carlo simulations to check our photom¬ 
etry at each wavelength: injection of sources in the data, and 
blind extraction, in order to measure the photometric accuracy 
(Fig. |3}. As expected, the photometric accuracy is of the order 
of 10% at flux densities larger than a few tens of mJy, and de¬ 
creases towards smaller flux densities closer to our noise level 


for nominal SPIRE map pixels of 6, 10, and 14", respectively. 
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250 /<m 350 yum 500/un 



^250 [mJy] 

^350 [mJy] 

^500 [mJy] 


Figure5. Differential Euclidean-normalized number counts, S 25 dN/dS , for various data sets (not corrected for flux-boosting or 
incompleteness) used to measure the relative behaviour. Symbols show: Planck IN region (red circles, 228 fields); Planck OUT 
regions (black circles); Planck whole field, i.e., IN+OUT (blue circles); sample of 535 z < 1 clusters of the HLS (open square); and 
HerMES Lockman SWIRE blank field (crosses). Our SPIRE sources corresponding to the Planck IN regions clearly show an excess 
in number density at 350 and 500 /mi, illustrating their red colour, and a potentially higher redshift than average. See Sect. 14. ll lor 
details. 


(dominated by confusion) at around 30 mJy. The completeness Table 3. 80 % completeness levels obtained from Monte Carlo 
levels are also measured, and reported in Fig.[4]and Table[3] source injection in the SPIRE maps for different fields, in mJy. 

Our Planck fields are denoted “ Planck high-z”. 


3.3. Band-merged catalogues 

We use the 350/rm map-based catalogue from StarFinder as 
the input catalogue. This wavelength has two advantages: first, 
its angular resolution meets our need to identify the sources; 
and second, it is consistent with the Planck colour selection, 
which avoids the many low-z galaxies peaking at 250/mi. As 
will be shown in Sect. 14.11 our sample has an excess in num¬ 
ber at 350/rm compared to reference samples drawn from the 
HerMES Lockman SWIRE level 5 field or the low-redshift HLS 
cluster fields, in agreement with the existence of an overdensity 
of SPIRE sources in the Planck beam. 

Our band-merging procedure has three steps. First of 
all, we optimize the measurement of the source position 
using the 250 //m channel, where available. Secondly, we 
measure a preliminary flux density on each source, which 
will serve as a prior to avoid unrealistic flux measurements 
while deblending. Thirdly, we perform spatially-simultaneous 
PSF-fitting and deblending at the best measured positions 
with the newly determined prior flux densities as inputs for 
FastPhot (iBethermin et al.l l2010bl) . The details are described 
in Appendix [A] This method provides better matched statis¬ 
tics (more than 90% identifications) than the blind extraction 
method performed independently at each wavelength (Sect. I3T21 
Table lAT2l vs. Table IaT1 >. We finally have a total of about 7100 
SPIRE sources, of which about 2200 are located in the Planck 
IN regions (giving an average of about 10 SPIRE sources per 
Planck IN field). 



80 % Completeness 

Field 

250 [t m 

350//m 

500 /rm 

Planck high-z. 

35.2 

37.0 

40.7 

HerMES Lockman-SWIRE_ 

35.0 

38.4 

Al.l 

HLS. 

49.3 

48.5 

54.3 


ing (Eddington bias), since we are interested in the relative be¬ 
haviour between the samples. The detected sources used here are 
extracted using the blind technique (Sect. I3T21) . We cut the sam¬ 
ples at 4 cr for these counts (<x values in Table [2j. We used five 
different data sets to estimate the number counts: 


Planck IN+OUT, our 228 Planck entire fields, each covering 
about 20' x 20'; 

Planck IN, our 228 Planck fields, using only the central parts 
corresponding to the Planck 50 % contour level (determined 
separately for each source), called the IN region; 

Planck OUT, our 228 Planck fields, using only the part exte¬ 
rior to the Planck 50 % level, called the OUT region; 
HerMES Lockman SWIRE, HerMES level"~5 field in 
Lockman ( Oliver et al.l[2010(). cov ering 18.2 deg 2 ; 

HLS, 535 cluster-fields of Egami et all (1201 Ol) . five from the 
Herschel Lens Survey KPOT, 282 from OT1, and 248 from 
OT2, hereafter refered to as HLS. 


4. Statistical Analysis of the Sample 

4. 1. Number counts: significant excess of red sources 

We compute the differential Euclidean-normalized number 
counts S 25 dN/dS , with S being the flux density at wavelength 
A, and N the number of sources per steradian. The counts 
are not corrected here for incompleteness, or for flux boost- 


The number counts of all these data sets are plotted in Fig. [5] 
and reported in Appendix [B] in Tables IB.II IB.2I and IB. 31 We 
derive from the total number count the following results. 

(i) The counts measured in the entire Planck fields (blue 
circles, meaning IN + OUT), HerMES Lockman, and 
HLS clusters, are generally compatible with each other at 
250 /mi, except for a few points at large flux density, above 
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300 mJy, where the numbers are small. This means that, 
on average, those fields do not show any strong deviations 
between them. 

(ii) The observed counts (uncorrected for incompleteness) 
show the characteristic shape of the Eddington bias: a cut¬ 
off below about 40 mJy, an excess around 50 mJy, and a 
behaviour compatible with the models at higher flux den¬ 
sities. It is beyond the scop e of this paper to re-derive un- 
biased number counts (e.g. , Glenn et alj|201 Ok [Oliver et al.l 
12010 b [Clements et 2010 ). as we only focus on the rela¬ 
tive trends. 

(iii) The IN counts at large flux densities (S > 300 mJy) show 
a systematic excess at 350 and 500 /am over the counts in 
the wide blank field (Lockman) as well as cluster fields. At 
300 mJy, the overdensity factors are, respectively, 4.1 and 
3.7 at 350gm, and about 20 and 16 at 500pm. This chro¬ 
matic excess (i.e., larger excess at longer wavelengths) is 
consistent with the hypothesis of the presence of a popula¬ 
tion of high-redshift lensed candidates at z = 2-A. We will 
show in Sect.[5]that this is indeed the case. 

(iv) The IN counts at 250pm are a bit higher (for S 250 < 
lOOmJy) than in Lockman, and are lower at larger flux 
densities. This means that the SPIRE counterparts of the 
Planck sources have, on average, counts that deviate little 
from blank fields at this wavelength. 

(v) The IN counts at 350 /am for S 350 — 50 mJy are a factor be¬ 
tween 1.9 and 3.4 higher than in Lockman and in the clus¬ 
ter fields. This means that, on average, the SPIRE coun¬ 
terparts of the Planck sources have a significant excess of 
350 pm sources compared to wide blank fields or z < 1 
HLS cluster fields. This is expected, given our Planck se¬ 
lection criteria. 

(vi) The IN counts at 500 gm for S 500 — 50 mJy are a factor of 
2.7-8 higher than in the Lockman and cluster fields. As in 
item (v) above, this means that, on average, the SPIRE im¬ 
ages of the Planck source targets show a significant excess 
of 500 pm sources compared to wide blank field or z < 1 
HLS cluster fields. 

(vii) The OUT counts are compatible with the wide blank extra- 
galactic fields, as well as the cluster fields. Our OUT zones 
can thus also be used as a proxy for the same statistics in 
blank fields. 

As a conclusion, the SPIRE observations of the Planck fields 
reveal “red” sources. The SPIRE images exhibit a significant ex¬ 
cess of 350 and 500 pm sources in number density compared 
with wide blank fields (HerMES Lockman SWIRE of the same 
depth) or fields targeting z < 1 galaxy clusters (HLS). This sig¬ 
nificant excess should be expected given the Planck colour se¬ 
lection, and is now demonstrated with secure SPIRE detections. 
It is therefore clear that there is no significant contamination by 
cirrus confusion in our Planck sample, and that where there is a 
Planck high-z candidate, Herschel detects galaxies. 


sources computed in the Planck OUT region, at SPIRE wave¬ 
length A. We have already shown (Sect. | 4 ~X 1 and Fig.0 that the 
OUT region has a density equivalent to that of blind surveys, and 
is thus a good estimate of p, the mean surface density. To reduce 
the Poisson noise, we use the counts from all the OUT regions. 

The overdensity contrasts b /( extend up to 10, 12, and 50 
at 250, 350, and 500 pm, respectively, with a median overden¬ 
sity 6 a of 6250 = 0.9, £350 = 2.1, and £500 = 5.0. This means 
that our Planck IN regions have an excess of SPIRE sources. 
Indeed, there are 50 fields with £500 >10, and 129 with bsoo > 4 
(with significance levels always higher than 4 cr^, see below). At 
350pm, there are 19 fields with £350 > 5, 37 fields with £350 > 4 
and 59 fields with A 350 > 3, as shown in Fig. [ 6 ] (left panel, in 
blue). In Appendix I dI we also use the densities measured with 
AKDE to estimate the overdensities. 

How significant are these overdensity contrasts? To quantify 
this we compute the mean density field using the AK DE a lgo¬ 
rithm (adaptati ve kernel density estimator , see Valtchanov et al.l 
120131 and also iPisanil 1 1 9961 : iFerdosi et al.l 1201 1 i. The principle 
is to generate a two-dimensional density field based on the po¬ 
sitions of the sources from a catalogue, filtered (smoothed) ac¬ 
cording to the source surface density. From this smoothed field, 
we compute the standard deviation and hence we derive the sig¬ 
nificance cr 6 of the overdensity. We also run 1000 Monte Carlo 
runs to get a better estimate of the scatter on cr 5 (by creating 
AKDE density maps using random source positions - but oth¬ 
erwise using the real catalogs of each field - and measuring 
the RMS over those 228000 realizations). All our fields show 
overdensities larger than 1.8 erg, with a median of 7 erg (Fig. [ 6 ] 
right panel). The typical number of sources in an overdensity is 
around 10 , but with a fairly wide scatter. 

We can then choose to select the reddest sources, this be¬ 
ing a possible signature of a higher redshift or a colder dust 
temperature. We define the SPIRE red sources with the follow¬ 
ing cuts in colour based on the distributions shown in Fig. [TJ 
S 350 / S 250 > 0.7; and S 500 /S 350 > 0.6. When selecting only the 
red sources, the overdensity significance increases: The mean is 
now 12 erg and the median 9 erg . We have 50 % of the sample at 
10 erg or more (when selecting the red sources), which is more 
than a factor of 3 larger than when selecting all SPIRE sources. 
23% of the sample is above 15cr,s (Fig. [6] right panel), i.e., 51 
fields. 

These high significance levels can be contrasted with the 
mean 6500 — 0.25 obtained by iRigbv et~akl (1201 4 ) at the loca¬ 
tions of 26 known protoclusters around very powerful radio¬ 
galaxies, drawn from the list of 178 radiogalaxies at z > 2 of 
iMilev & de Breuckl (12008 ). They can also be c ompared to the 
few similar examples in Clements et alJ (12014) . with at most 
4.7 erg at 350 pm. We show in Fig. [6] right panel, in orange, the 
cumulative normalized significance erg of 500 random positions 
in the Lockman field; Using this test sample illustrates the high 
significance of the overdensities of our sample. 

Examples of some overdensities are shown in Appendix [E] 
where we present a gallery of representative SPIRE data. 


4.2. Overdensities 

We compute the dimensionless overdensity contrast 62 of our 
fields at wavelength A via 

„ Pin-Pout ... 

oa = -, ( 1 ) 

Pout 

where /pxj is the surface density of SPIRE sources in the Planck 
IN region, and pout is the mean surface density of SPIRE 


4.3. Colours of the sources 

Using the band-merged catalogue (Sect. [Oi l, we can derive the 
colours of the sources, i.e., the ratios of the observed flux densi¬ 
ties S500/S350 vs. S 250/S350- In this colour-colourspace, redder 
sources will have higher S 500/S 350 ratios and lower S 250 /S 350 
ratios, and will tend to lie in the lower left part of the diagram. 

We provide in Fig. [ 7 ] the source surface density histograms 
for the two SPIRE colours: S 350/S 250 and S 500/S 350 (red, Planck 
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Figure 6. Left : Cumulative histogram of the overdensity contrast d’350 (blue) of each Planck field (based on 350 /am SPIRE sources). 
Overdensities are fairly large, with 59 fields having £350 > 3 . Right: Cumulative (normalized) statistical significance in cr derived 
from the density maps. Blue represents all our SPIRE sources, red represents only redder SPIRE sources, defined by S 350/.S' 250 > 0.7 
and S 500/ S 350 > 0 . 6 , and orange 500 random fields in Lockman. Most of our fields have a significance greater than 4 <x, and the 
significance is higher still for the redder sources. See Sect. l 4 . 2 l for details. 


to 

c 

a> 

"O 

CD 

o 
1 — 
3 

o 

CD 



0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 


2.5 3.0 


^ 350^250 


^ 500^350 


Figure 7. Colour counts: source surface density as a function of the SPIRE colour, S 350 /S' 250 (left), and S 500 /S 350 (right). Histograms 
are: red solid line, Planck IN; black line, Planck OUT; green, Planck IN lensed fields only; blue line, z < 1 HLS clusters; and 
orange dashes, Lockman SWIRE. The Planck IN sources (total and/or lensed sources) show a much higher surface density than 
other samples, owing mainly to our all-sky search strategy. See Sect. l 4 . 3 l for details. 


IN; black, Planck OUT; green, only the lensed fields; blue, HLS; 
and dashes, HerMES). The IN sources show three times larger 
surface densites in S 350 /S 250, and four times larger surface den¬ 
sities in S 500/.S' 350 than OUT and Lockman sources. The IN 
source distribution peaks at much higher surface density than 
any other sample, suggesting that our sample is dominated by 

red and overdense SPIRE sour ces. _ 

Following the approach of lAmblard et ah ! ( 2010 ). we gener¬ 
ate the SED of 10 6 modified blackbodies with 10 % Gaussian 
noise and explore three parameters: the blackbody temperature 
( T) in the range 10-60 K; the emissivity (J 3 ) in the range 0 - 2 ; 


and the redshift for z = 0 - 5 . For each set of parameters we 
convert the fixed, observed SPIRE wavelengths into rest frame 
wavelength (at redshift z, thus varying with z) using: T re st = 
4spire/(1 + z). We then calculate the flux at each wavelength 
using T rest and compute the colours S 250 /S 350 and S 500/S 350 for 
each set of parameters. Fig. [8] shows those realizations, colour- 
coded in redshift (deep purple for z = 0 , up to red for z = 4 ). 
Our SPIRE IN sources are shown in black. They predominantly 
fall in the redshift region corresponding to z ~ 1 . 5 - 3 . We note, 
however, that the redshift-temperature degeneracy is present in 
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Figure 8. Herschel-SPVRE colour-colour diagram of IN sources (black dots) on top of redshifted blackbodies. Left : colour codes 
redshift. Right: colour codes dust temperature in kelvin. Although there is a degeneracy between redshift and temperature, the 
colours suggest a range of z ~ 1.5-3 for the SPIRE sources corresponding to the Planck IN regions. The upper right symbol gives 
the typical (median) error bars of the measured points. See Sect. l4.3l for details. 



Maximum flut density tmjy) 


Figure 9. Classification of 228 Planck fields with Herschel- 
SPIRE. The x-axis represents the maximum flux of the SPIRE 
source within the IN region: lensed source candidates are se¬ 
lected if they fall above 400 mJy (i.e., the right column, seven 
sources). The v-axis represents the wavelength at which the 
brightest SPIRE source peaks in the IN region. Overdensities are 
thus selected in the lower left four cells. The colour represents 
the number of fields in each cell (as shown in the colour bar), 
the numbers in the cells are the number of fields and percent¬ 
age in each cell. Our sample is thus dominated by overdensities 
peaking at 350/rm. See Sect. l4.4l for details. 


this suggest ive result (lAmblard et al.l20TotlPope & Charvl2010t 
iGreve et aJ1l20 1 2l) . and will be discussed in Sect. 16.31 



Figure 10. A high-z cluster candidate observed by Planck, 
Herschel, and Spitzer- IRAC (image covering 0.5' x 0.3'). We 
show the IRAC channel 1 (3.6 pm) image, with SPIRE 350/rm 
white contours overlaid. Colour contours represent statistical 
significance of the local overdensity, from light blue to red: 3, 
4, 5, 6, and 7 erg. See Sect. l6.1l for details. 


4.4. Classification of the sources 

We classify the Herschel fields in an automated way and for 
preliminary analysis into two main categories: (1) overdensities, 
appearing as clustered sources, and (2) lensed candidates, ap¬ 
pearing as a bright single source as a counterpart of the Planck 
source. To make this classification, we measure within the IN 
Planck region the flux density and colour of the brightest SPIRE 
source. 
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We sort the fields using two criteria: (1) the SPIRE chan¬ 
nel at which the source has its maximum flux density; and ( 2 ) 
its flux density at 350 /am. We will use this two-parameter space 
to classify the sources. This classification provides the first set 
of information about the sources, although a definitive classifi¬ 
cation will require follow-up data to confirm the nature of the 
sources. 

In this space, the lensed candidates will populate the “red” 
and “bright” areas (typically 350 or 500gm peakers, and S v > 
400 mJy). This does not necessarily mean that all the sources in 
this area will be strongly gravitationally lensed (although this is 
confirmed by further follow-up, see following sections). 

The overdensity candidates will be located in a different area 
of this two-parameter space than the lensed candidates. They 
will tend to populate the fainter end (typically S v < 200 mJy) 
and will have red colours (typically 350 or 500gm peakers). 

We do not expect many 250 pm peakers because of the 
Planck colour selection, for which there is a bias towards red¬ 
der colours and thus potentially higher-redshift sources. 

Fig. [9] summarizes the classes of SPIRE identifications of 
the Planck high-z (meaning z > 1.5) candidates, based on those 
criteria. We show that the vast majority of our sample is com¬ 
posed of overdensities of 350pm peakers. The number of lensed 
candidates is rather small in comparison, amounting to seven 
sources based on this criterion. All are confirmed to be strongly 
gravitationally lensed galaxies (ICanamerasl2015l) . We will show 
later that visual inspection led us to discover a few more lensed 
sources. 

The dashed lines in Fig. |7| show the colour vs. surface den¬ 
sity of the lensed candidate population, whose red colours are 
confirmed. 

Finally, we derived number counts (as in Fig. 0 but now by 
separating the overdensities from the lensed sources that dom¬ 
inate the counts at large flux densities, typically S > 250 mJy 
(shown in Appendix [C] in Fig. 1C.II ). This separation suggests 
that the excess at large flux densities in the number counts of our 
sample is due to the presence of bright lensed sources, compared 
to reference samples (HLS, HerMES). 


5. Strongly gravitationally lensed source 
candidates 

5.1. Validation of existing lensed sources 

INegrello et al.l (120071) and iBethermin et al.1 (1201 2l) predicted 
that a small, but significant fraction of very bright high- 
redshift (z > 2 ) submillimetre galaxies are strongly grav¬ 
itationally lensed, dusty starbursts in the early Universe. 
INegrello et al.l (1201 Ob pres ented obs e rvatio nal evidence of these 
predictions. In addition, iFu et al.l (12012b confirmed the na¬ 
ture of the source H-ATFAS J114637.9-001132 as a strongly 
gravitationally lensed galaxy at z — 3.3; this source was 
part of the first release of th e Planck ERCSC Catalogu e 
(IPlanck Collaboration VIIl 1201 It iPlanck Collaboration! 1201 ll) . 
and fell fortui tous ly into the H-ATFAS survey field (see also 
iHerranz et al l l201 3». Another independent confirmation that 
Planck sources can be strongly gravitationally lensed came 
from the source HLS J091828.6+514223, confirmed as a bright, 
z = 5.2 gravitationally lensed galaxy behind the massiv e 
intermediate-redshift galaxy cluster Abell 773 (ICombes et al. 
2012), as part of the Herschel Lensing Survey (lEgarni et al. 
20101) . This source was independently found in our survey and 
is part of the Herschel /SPIRE OT2 selection. It is excluded from 
our sample to satisfy the condition of non-redundancy. 


5.2. Previously unknown gravitationally lensed sources 


In the absence of extensive follow-up, it is challenging to distin¬ 
guish between single or multiple strongly gravitationally lensed 
galaxies behind the same foreground structure and overdensi¬ 
ties of intrinsically bright submillimetre galaxies, for all but the 
brightest gravitationally lensed sources. Moreover, at flux den¬ 
sities of about S 350 = 250 mJy and above, isolated SPIRE point 
sources may turn out to be associations of multiple FIR g alax ies 
when observed at higher spatial resolution (Ivison et al.112013l) . 

As a first step towards identifying the gravitational lensed 
candidates in our sample (Sect. 14.4b . we therefore focused on 
those targets where SPIRE shows only a single, very bright 
source, with the typical FIR colours of high-redshift (z > 2) 
galaxies associated within one Planck beam. We thus identified 
seven isolated SPIRE point sources, as described in Sect. 14.41 
plus five others at slightly fainter flux densities. All have peak 
flux densi ties at 35 0 um, i ncluding 11 with S 350 = 300- 
1120 mJy (ICanamerasll20i5 ). Six of these galaxies were taken 
fr om the PCCS, the remainin g six originate from the sample 
of lPlanck Collaboration! ( 120151) . Although the initial selection of 
lensed candidates was carried out by eye upon the reception of 
the SPIRE imaging, seven of these targets were also identified 
with our automatic classification (Sect. l4~4l >. We used the IRAM 
30-m telescope to obtain firm spectroscopic redshifts via a blind 
CO line survey with the wide-band receiver EMIR. We iden¬ 
tify 2-6 lines per source, which confirms they are at redshifts 
z — 2.2-3. 6 . Interferometry obtained with the IRAM Plateau 
de Bure interferometer and the Submillimeter Array, as well 
as several empirical calibration s based on FIR luminosity and 
dust temperature (following lHarris et+if 2012^, and CO line lu¬ 
minosity and line width (following iGreve et alj|2012h . demon¬ 
strate that these are indeed strongly gravitationally lensed galax¬ 
ies, amongst the brightest on the submm sky ( ICanamerasll2015l) . 


6. Candidate high-z overdensities 

Without rejecting the possibility that the observed overdensities 
of red SPIRE sources could be chance alignments of structures 
giving coherent colours (e.g. JChiang et al.|l2013l) . we can pon¬ 
der the nature of those overdensities. Could they be high-redshift 
intensively star- forming galaxy pro t o-clusters? I n deed, recent 
studies, such as Gobatet al. J 201 ll) . ISantos et al.l (1201 ll 1201 3l 
120141) . and IClements et al.l (1201 4l) confirm the presence of high 
redshift (z > 1.5) galaxy clusters emitting enough energy in the 
submillimetre to be detected. Those previously detected clusters 
are, however, in a different state of evolution than our candidates. 
Many of the confirmed clusters exhibit X-ray emission, suggest¬ 
ing already mature and massive clusters. The following sections 
investigate different aspects of our sample, which is uniquely se¬ 
lected by strong submm emission over the whole sky, enabling 
us to unveil a rare population. 


6.1. First confirmations 

It is beyond the scope of this paper to summarize all the follow¬ 
up observations conducted so far, but we give here three high¬ 
lights, which will be discussed in more detail in subsequent pa¬ 
pers. 

Our sample contains a source from Herschel OT1 that was 
followed up early as a pilot programme. Using CFHT and Spitzer 
imaging data, and VLT and Keck spectroscopy, a structure was 
identified at z ~ 1.7-2.0 (Flores-Cacho et al., 2014, in prep.). By 
fitting a modified blackbody and fixing redshifts at respectively 
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Figure 11. Left. SPIRE 3-colour image of one Planck high-z overdensity candidate (located in the centre), in an image covering 
about 20' x 20'. We note the bluer source on the left (circled), which cor respon ds to the galaxy cluster shown on the right. Right: 
a confirmed z = 1-58 galaxy cluster XMMU J0044.0-2033 (Santos et al.|[201 lh which lies close to our Planck high-z source. The 
background image is a two colour composite of Spitzer -IRAC channel 1 (3.6 pm) and channel 2 (4.5 pm), covering about 3' x 2'. 
The contours are Herschel- SPIRE 350 pm. See Sect. 16. 1 1 for details. 



(b) HLS sample 


1 trcmtn 


(d) Random in HarMES 


(c) HarMES sample 


Figure 12. Stacks of SPIRE 350/rm data (8.7' x 8.7'): (a) 220 Planck Herschel fields; (b) 278 HLS clusters; (c) 500 sources in the 
HerMES Lockman field, peaking at 350 pm: and (d) 500 random positions in the HerMES Lockman field. Black contours show the 
density of red sources (using AKDE). Our Planck fields clearly exhibit a high significance of extended (a few arcminutes) submm 
emission, due to the presence of many red point sources, not seen in the HLS or HerMES sources. See Sect. l6.2l for details. 
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1.7 and 2.0, we find an average dust temperature for the sources 
of, respectively, 27K and 35K. 



Redshift [z] 


Lockman field; and (d) 500 random positions in the HerMES 
Lockman field as a null test. The stacks are presented in Fig.fl2l 


We can see that the Planck fields show clear and signifi¬ 
cant extended (a few arcminutes) emission due to the cluster¬ 
ing of bright submm sources. This kind of extended emission 
of clustered submm sources is not observed in the HLS clus¬ 
ters, nor around 350pm peaking HerMES sources. The ran¬ 
dom stack_vafidatesthe_absence_ofa_systematic_effecthi stack¬ 
ing ( Dole et aflbood Heth er min et al.ll2010at lYiero et aPl2013b 


iPlanck Collaboration XVIIlll2014l) . This proves that our sources 

are of a different nature than mature HLS clusters or average 
HerMES submillimetre sources. 


We overplot in Fig.[l2]the contours in density (see Sect. 14.21) 
of the SPIRE red sources (defined in Sect. 14.21 ). The stacks indi¬ 
cate that the Planck sample exhibits strong overdensities of red 
SPIRE sources. By contrast, the HLS sample shows a smooth 
and weak density of red sources (except for the presence of some 
point sources - the background lensed galaxies). As expected, 
the HerMES sample shows a stack consistent with the SPIRE 
PSF, while the random sample is consistent with noise. 


Figure 13. SPIRE photometric redshift distribution of the 
roughly 2200 SPIRE sources in the IN regions, as a function 
of the fixed assumption for dust temperature: 7 d = 25, 30, 35, 
40, and 45 K (from left to right). See Sect. l6.3l for details. 


Another field has been imaged by Spitzer as part of our GO- 
9 programme - see Fig.|T0]for the IRAC 3.6 pm image with the 
SPIRE 350 pm contour. The colour ratios of IRAC flux densi¬ 
ties 3.6pm/4.5 pm of the sources detected in the overdensity 
exhibk red_colours, indicating their probable high-redshift na¬ 
ture (lPapovichll2008l) . The overdensity of IRAC 3.6pm sources 
has a statistical significance of 7 cr s (Martinache et ah, 2014, 
in prep.). While not yet confirmed as a proto-cluster or cluster 
through spectroscopy, this structure, seen with Spitzer at near- 
infrared wavelengths, supports the hypothesis of strongly clus¬ 
tered highly star-forming infrared sources. 

Finally, we have also found one SPIRE field with overlap¬ 
ping d ata from the literature. The observations from lSantos et al.l 
(1201 ll) confirm the presence of a z = 1.58 cluster using op- 
tical/NIR spectroscopy. This particular high-z cluster has also 
been detected in X-rays, leading to a total mass estimate of about 
3-5 x 10 14 M o . It is located a few arcminutes from one of our 
targets (Fig. |TT]right), so is not formally an identification of our 
Planck source, but rather a coincidence. While not connected to 
our target source here, this cluster can nevertheless bring use¬ 
ful information. The SPIRE colours of this cluster are bluer than 
our sources. Since our target shows redder SPIRE colours, this 
suggests it has a higher redshift or a cooler dust temperature. 


6.2. Large clustering of our sample revealed with stacking 

We can investigate if the observed overdensities have differ¬ 
ent clustering properties than two test samples: (1) the HLS 
massive z < 1 clusters; and (2) 350pm peaker sources in 
the HerMES Lockman field brighter than 50mJy. To do so, 
we use a stacking analysis (see, e. g., iMontier & Giardl 1 20051 
|Dole_eta2 2006^ Bra glia et al.ll201 ll) . here applying the method 
of iBethermin et al.l ( 2010a ). We stack the following SPIRE 
350 pm data: (a) the 220 overdensities of our sample; (b) 278 
HLS clusters; (c) 500 bright 350pm peaker sources in the 


6.3. Dust temperatures, photometric redshifts, luminosities, 
and star formation rates 


With the hypothesis that our overdensities are actually gravi¬ 
tationally bound structures, i.e., clusters of star-forming, dusty 
galaxies, the overdensity colours (Fig. [7]) suggest a peak redshift 
around z = 1.5-3 (Fig. [8])- Using a modified blackbody fit (with 
/? = 1.5) and fixing the dust temperatures to 7j = 25, 30, 35, 40, 
and 45 K, we obtain the redshift distributions shown in Fig. [13] 
for the roughly 2200 sources found in the IN regions for these 
overdensities. For T& = 35 K, the distribution peaks at z — 2. For 
higher dust temperatures, the peak of the distribution is shifted 
towards higher redshifts, since dust temperature and redshift are 
degenerate for a modified blackbody, where 7^/(1 +z) =constant. 


Greve et al 


Many studies (e.g.. iMagdis et al' 1201 Ot Elbaz et al 


Weiss et al. 


20421 Magdis et al . 2012; Svmeonidis et al 

if 


2011 


2013 


2013t iMagnelli et al.l 20141) suggest that dust tern 


peratures of z ~ 2 sources are typically of the order of 35 K, 
consistent with the few measurements in hand (see Sect. I6.lt 
With the conserva tive as sumption that T& = 35 K for all sources 
(e.g., iGreve et al.l 120121) . we can derive the IR luminosity for 
each SPIRE source (by integrating between 8 pm and 1 mm), 
and find that it has a broad distribution peaking at 4 x 10 12 L G 
(Fig. IT4] i. Assuming that the conventionally a ssume d relation¬ 
ship between IR luminosity and star formation (iKennicutdl 1998b 
lBellll2003h holds, and that AGN do not domi nate (as was also 
fou nd for the objects prev iously discussed bv ISantos et al.ll2014l 
and lClements et akll2014l) . this would translate into a peak SFR 
of 7OOM 0: yr 1 per SPIRE source (Fig. [T5j. If we were to 
use colder dust temperatures, the peak SFR would be at about 
200 M 0 yr 1 per SPIRE source (Fig. [T31 . If we were to use 
warmer dust temperatures, the SFR would increase because the 
implied redshift would be higher and still compatible with other 
o bservati ons, e.g., the compilation of submillimetre galaxies in 
iGreve et akl (12012l) . where Lir is measured in the range 3 x 10 12 
1O 14 L 0 and T& in the range 30-100K, for redshifts above 2. 


We can also estimate the IR luminosities and SFRs of the 
overdensities, i.e., potential clusters of dusty galaxies at high-z, 
rather than the single SPIRE sources that we used above. With 
typically 10 SPIRE sources in the IN region, we sum up the IR 
luminosities of SPIRE sources within the IN region (under the 
same assumption of - 35 K) to obtain a peak IR luminosity 
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of 4 x 10 13 L 0 . This translates to a peak SFR of 7 x 10 3 M 0 yr -1 
per structure, as seen in Fig. [16] 

Our estimates of Lir are within the range of FIR luminosi¬ 
ties expected fo r massive, vig orously star-forming high-z struc¬ 
tures (e.g.. iBrodwin et afll2013l) . which are perhaps protoclus¬ 
ters in their intense star formation phase, and are c onsistent 
with the four bound structures found in lClements et al.l (120141) . 
These z > 2 intensively star forming proto-clusters are also ex¬ 
pected in some models, e.g., the proto-spheroids of ICai et~akl 
(1201 3l). or the halo s harbo uring intense star formation discussed 


in 


Bethermin et al.l (120131) . 



IR Luminosity per source [L 0 ] lel3 

Figure 14. Infrared luminosities (Lir) of the approximately 2200 
SPIRE sources in the IN regions, as a function of the fixed dust 
temperature: 35 K (purple, right histogram), and also shown for 
illustration 25 K (blue, left histogram). See Sect. l6.3l for details. 


6.4. Massive galaxy clusters in formation? 

It is of course tempting to postulate that all (or at least a subset) 
of our sources include massive galaxy clusters in the process of 
formation. However, demonstrating this conclusively with the 
present Herschel photometry alone is not possible. We have 
therefore undertaken a comprehensive multi-wavelength photo¬ 
metric and spectroscopic follow-up campaign to directly con¬ 
strain the nature of our sources and in particular, to investigate 
whether they are good candidates for being the progenitors of to¬ 
day’s massive galaxy clusters seen during their most rapid phase 
of baryon cooling and star formation. 

Even without such an explicit analysis of the astrophysical 
nature of our candidates, we can already state that: ( 1 ) our esti¬ 
mated luminosities are consistent with the expected star forma¬ 
tion properties of massive galaxy clusters, as obtained from stel¬ 
lar archaeology in galaxies in nearby clusters and out to redshifts 
Z ~ 1; and (2) our total number of overdensities are consistent 
with the expected numbers of massive high-z galaxy clusters, as 
obtained from models of dark matter halo structure formation. 

For point (1), a broad consensus has now been developed 
whereby most of the stellar mass in massive cluster galaxies 
was already in place by z — 1. The most stringent observational 
constraint is perhaps the tight red sequence of cluster galaxies 
in colour-magnitude diagrams of massiv e galaxy clusters in the 
optical and near-infrared (e.g.. IRenzinil 12006 ). which suggests 


that massive cluster galaxies formed most of their stellar popula¬ 
tions within a short timescale of < 1 Gyr, i.e., within one cluster 
dynamical time, with little star formation thereafter. Rest-frame 
optical studies of distant clusters suggest that the bright end of 
the red sequence was already in place by z - 1 (iRudnick et al.l 
120091) . with a likely onset amongst the m ost massive galaxies by 
- 2, or even before dKodama et al.1120071) . 



SFR per source [M Q yr 1 ] 


Figure 15. Star formation rate (SFR) for individual SPIRE 
sources in the overdensity sample (i.e., in the Planck IN region). 
For 7’,| = 35 K, a dust temperature favoured by the literature, 
(purple, right histogram) the distribution peaks at 7OOM 0 yr~ 1 . 
We also show for illustration 7)| = 25 K (blue left histogram 
peaking at 2OOM 0 yr '). See Sect. l6.3l for details. 


The shortness of the star formation period in nearby clus¬ 
ters allows us to derive an order-of-magnitude estimate of the 
total star formation rates during this period, as implied by the 
fossil constraints, and the resulting FIR fluxes. To obtain upper 
limits on the stellar mass formed during this epoch, we use stel¬ 
lar mass estimates obtained for 93 massive X-ray-selected clus¬ 
ters w ith M 500 = (10 14 -2 X 1O 15 )M 0 at z = 0-0.6 bv lLin et al.l 
(1201 2l) based on WISE 3.4/im photometry. M 500 is the total mass 
within the central cluster regions where the mass surface den¬ 
sity exceeds the cosmological value by at least a factor of 500. 
iLin et al.l find that in this mass range, stellar mass scales with 
M 500 as (M*/10 12 M o ) = (1.8 + 0.1) (M 5 oo/lO 14 M 0 )°- 71±ao4 . For 
clusters with M 500 = 1 X 1O 14 M 0 (or 2 X 1O 15 M 0 , their highest 
mass), this corresponds toM, = 2 x 1O 12 M 0 (or 1.5 x 1O 13 M 0 ). 
We do not expect that including the intracluster light would in¬ 
crease th ese estimate s by more than a few tens of percent (e.g., 
iGonzalez et al.ll2007t) . which is a minor part of the uncertainty in 
our rough order-of-magnitude estimate. 

These mass and timescale estimates suggest total star forma¬ 
tion rates in the progenitors of massive galaxy clusters at lower 
redshifts of a few x 10 3 and up to about 2xl0 4 M 0 yr 1 . Using the 
conversion of lKennicutti ( 1998 1 between star formation rate and 
infrared luminosity (Sect. ROK this corresponds to Lir - (10 13 - 
10 14 )L o . This is well within the range of the global star forma¬ 
tion rates and Lrir values in our sample of overdensities. 

For point (2), the absence of such objects, we use the 
[Tinker et al.l (120081) halo model to compute the expected sur¬ 
face density of dark matter halos. We expect between 8500 and 
1 x 107 dark matter halos at z > 2 having masses M tot > 10 14 M 0 
and 10 13 M o , repectively, over 35% of the sky. With the as- 
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Total SFR [M 0 yr '] 


Figure 16. Total star formation rate (SFR) per Planck IN region, 
i.e., for the overdensity sample, assuming each overdensity is 
an actual high-z cluster of dusty galaxies, located at a redshift 
given by the SPIRE photometric redshift, and using a modified 
blackbody fit with a fixed dust temperature of T& — 35 K (purple, 
right histogram peaking at 7 x 10 3 M G yr -1 ), this temperature 
being favoured by the literature. We also show the total SFR 
for T c i = 25 K, giving a peak at 2 X 10 3 M G yr 1 . The structures 
could be intensively star-forming clusters. See Sects. I631 and l6. 41 
for details. 


sumption that the cluster galaxies have formed most of their 
stellar populations within a short timescale (see above), we esti¬ 
mate that only a fraction of those halos will be observationally 
caught during their intense star formation phase, when they are 
infrared- and submm-bright. Assuming that this phase happened 
mainly between z - 5 and z = 1.5-2.0, we find that the upper 
limit on the formation period of each individual cluster, about 
1 Gyr, is 3-4 times shorter than the cosmic time elapsed over this 
epoch, 2-3 Gyr. We would therefore expect to find a few thou¬ 
sand dark matter halos on the submm sky at any given observing 
epoch. This order-of-magnitude estimate is above with our find¬ 
ing of about 200 overdensities with Planck in this study, and 
with a few hundred to be de tected in the full Planck data set by 
IPlanck Collaboration (201 5). but not more than an order of mag¬ 
nitude, illustrating the need to understand the detailed processes 
of star formation in the most massive halos (iBethermin et al.l 
1201 31 ) . 


7. Conclusions 

Our Planck sample based on a colour selection of cold sources of 
the CIB is overwhelmingly dominated in number by significant 
galaxy overdensities peaking at 350 //m, and a minority of rare, 
bright z > 1.5 strongly gravitationally lensed sources, among 
which are the brightest ever detected in the submillimetre. This 
confirms the efficiency of Planck to select extreme and rare 
submillimetre sources over the whole sky, as well as the need 
for higher angular resolution imaging using Herschel (and cur¬ 
rent/planned ground-based submm/mm observatories) to iden¬ 
tify and study them. 

From the analysis of Herschel observations of the sample of 
228 Planck cold sources of the CIB, we draw several conclu¬ 
sions. 


- Less than 2 % of the fields are Galactic cirrus structures. 
When a clean and controlled selection is performed on 
Planck , there is thus no reason to expect a large cirrus con¬ 
tamination when using a conservative Galactic mask. 

- With about 93 % of the overdensities peaking at 350 /am, and 
3.5 % peaking at 500yum, our sample unambiguously selects 
sky areas with the largest concentrations of red SPIRE galax¬ 
ies, consisting of highly significant overdensities. 

- Some of the overdensities are confirmed high-z structures, 
e.g., a source at z — 1.7-2.0 (Flores-Cacho et al. 2014, in 
prep.). 

- The significance in density contrast of the overdensities (e.g., 
half of the fields above 10cr 5 when selecting red SPIRE 
sources) is higher than any other sample targeting proto¬ 
clusters or high -z clusters in the submillimetre, confirming 
the relevance of the strategy of using Planck data on the 
cleanest parts of the sky to uncover high-z candidates. 

- The SPIRE sources in the overdensity fields have a peak 
redshift of z = 2 or 1.3, if we fix the dust temperature at 
7’,i = 35 K or 25 K, respectively. 

- With the assumption of T& — 35 K, each SPIRE source has 
an average IR luminosity of 4 x 10 12 L e , leading to star for¬ 
mation rates for each source peaking at 700 M G yr 1 . If con¬ 
firmed, these exceptional structures harbouring vigorous star 
formation could be proto-clusters in their starburst phase. 

- Assuming the SPIRE sources are located in the same 
large-scale overdensity, we derive a total IR luminosity of 
4 x 10 13 L o , leading to total star formation rates of 7 x 
10 3 Mr., yr -1 , and with around 10 detected sources per struc¬ 
ture. 

- About 3 % of our sample is composed of intensely gravita¬ 
tionally lensed galaxies. This sample is unique, as it targets 
the brightest observed such sources, typically above 400 mJy 
at 350//m and reaching up to the jansky level. They are all 
spectroscopically confirmed to lie at redshifts z = 2.2-3.6 
(ICanamerasll2015b lN esva dbafl20T5h . 

- The novelty and efficiency of our new sample is that it 
provides about 50 times more fields for a Planck-Herschel 
co-analysis than in existing Hers chel su rve ys searc hing for 
serendipitous Planck sources (e.g JClements et alJ 2014 ), for 
only a fraction of the Herschel observing time. 

- Our new sample exhibits high density contrasts with a high 
significance: for the red sources, 30 % of our sample (about 
70 fields) shows significance levels higher than 4.5 cr a , and 
15 % of our sample (34 fields) are higher than 7 cr 6 . 


The presence of the overdensities (and the hypothesis that 
these could be actua l forming, bound s t ructures) is consistent 
with other findings ( Santos_etl^_ 2014 JClements et al.l 12014b 


with other findin gs (laantos et alJ zUI4t lllements et al.l l/U14t 

iNoble et aD 12013b iWvlezalek et al. 12013bh . Our Planck extrac¬ 


tion and confirmation with Herschel provides an efficient selec¬ 
tion over the whole sky, biased towards star-forming sources, 
giving a few hundred candidates (Fig. Ell. This selection 
nicely complements the structures found at z ~ 2 already de¬ 
tected by different means, e.g., with the X-ray signature of 


in association with radiogalaxies (Pentericci et al. 1993, 

2000; 

Brodwin et alJ 2005: 

Milev et al. 2006; Nesvadba et al. 

2006; 

Venemans et al. 2007 

Dohertv et al. 201C|Pi 

inovich et al. 201C; 

Brodwin et afJ^OlcT 

201 lb Hatch etalloil 

Gobat et al. 

2011, 

Santos et al. 2011; 

Stanford et al. 201^ 

Santos et al. 

2013 

20 ll 

Brodwin et al. 

2013; Galametz et al. 

2013; Rigbvetal. 

20 ll 

Wylezalek et al. 2013a; Chiang et al. 

20 ll Cooke et al. 

2014b Cucciati et al. 2014). However, it is hard to estimate the 
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fraction that may b e comp rised of random alignments of unre¬ 
lated clumps (e.g. JChiang et al.ll2013l) . since the fluctuation field 
due to clustering is highly non-Gaussian. Ancillary data (partic¬ 
ularly redshifts) are needed to confirm the associations, and an 
ambitious follow-up programme is underway. 

To further characterize the overdensities discovered by 
Planck and confirmed by Herschel , we also started to stack them 
in many ancillary data sets . The first data set is the Planck lens- 
ing map (Planck Collaboration XVIIl2014l) . The goal is to detect 
a signature of the presence of large gravitational potential, as has 
been demonstrated already by correlating theJTIB fluctuations 
with the lensing map dPlanck Collaboration XVHll2014l) . The 
second data set is to use the nine frequenc y maps and the diffuse 
tSZ map obtained by MILCA dHurier et al.l2013l) . The goal here 
is to det ect a signatu re of the thermal Sunyaev-Zeldovich (tSZ) 
effect dSunvaev & Zeldovichll 19691. Il972lf and thus the presence 
of hot gas. These analyses will be published in a subsequent pa¬ 
per. 

The Planck data thus provide unique and powerful samples 
to uncover rare populations, exploitable for extragalactic studies, 
as was shown with our successful follow-up with Herschel of 
228 Planck high-z candidates. There are hundreds more Planck 
high-z candidate s (the full catalogue is bein g finalized and will 
be published in iPlanck Collaboration! 12015b to be studied and 
characterized. 
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Appendix A: SPIRE band-merging procedure 


Table A.l. Percentages of SPIRE sources matched at 250 and 
500 pm from the 350 pm input SPIRE catalogue, and the fre¬ 
quency of matches (i.e., number of sources found per 350/rm 
source). Here, IN refers to the Planck source region, and 
OUT refers to the zone outside the Planck 50% contour (see 
Sect. 12.3b . The rows are: no match (top); one match; two 
matches; and three matches (bottom). More than 60% of the 
SPIRE 350 p m sources corresponding to the Planck sources (i.e., 
the IN regions) have counterparts at 250 and 500 /mi. 


250 pm 500 gm 

IN OUT IN OUT 


No match . 10.9 16.7 38.8 67.0 

1 match. 82.0 77.0 60.7 32.9 

2 matches. 6.9 6.0 0.5 0.1 

3 matches. 0.1 0.2 0 0 


We describe here the details of the band-merging procedure, 
highlighted in Sect. 13.31 

First step: positional optimization of blind catalogues using 
shorter wavelength data. 

- Apply a 4 cr cut to the 350 pm blind catalogue obtained with 
StarFinder (Sect. l3T2l >. This cut is imposed to obtain a high 
purity of detections. 

- Match by position the 350 pm sources to the 250 pm cata¬ 
logue, within a radius of one 250 pm FWHM. 

The following depends on whether there are 0, 1, 2, or more 
matches: 

- if zero sources are found at 250 pm, this is a non-detection at 
250/tm; 


17 











Planck Collaboration: high-z sources discovered by Planck and confirmed by Herschel 


Table A.2. After the SPIRE band-merging process, the percent¬ 
age of SPIRE sources that do not have any flux density measure¬ 
ment (i.e., are blended or not detected or have multiple matches 
in Table IA. 11 1 . More than about 90% of the Planck IN sources 
have counterparts at all three SPIRE wavelengths. 



250 pm 

350 pm 

500 pm 


IN 

OUT 

IN 

OUT 

IN 

OUT 

No detections. . . 
S/N < 3 . 

. 0.2% 

. 3.2% 

0.6% 

7.6% 

0.2% 

0.9% 

0.6% 

3.2% 

0.5% 

10.6% 

3.6% 

39.4% 


- if one source is found, we take this source to be the 250//m 
counterpart; 

- if two sources are found, we replace the position of the 
350 pm unique source by the two positions of the 250 pm 
sources (the 350 gm flux density will be measured later) and 
take the two 250/mi flux densities; 

- we apply the same procedure for triply or more matched 
sources: we keep the positions of the 250 pm sources as pri¬ 
ors (only one case, outside the Planck beam. 

We repeat this operation at 500 gm using one 350 /mi FWHM as 
the search radius. The statistics of those basic positional matches 
are reported in Table IA.1I We have more than 60 % unique 
matches in the IN regions, i.e., for the Planck sources. This sim¬ 
ple method is, however, biasing the sample towards strong and 
unique detections at the three SPIRE wavelengths. The two fol¬ 
lowing steps address this issue. 

Second step: prior flux density determination. 

- At 350 pm: take the StarFinder flux densities (Sect. l3T2l > 
for single-matched source with 250 pm; assign S 350 = 0 if 
there are two or more matched sources at 250 pm. 

- At 250 pm: take the flux density for single-, double- or triple- 
matched sources with the 350 pm source; assign S 250 - 0 if 
there is no match with 350 pm at this step. 

- At 500 pm: take the StarFinder flux densities (Sect. l3T2l > 
for single-matched source with the 350 pm source; Assign 
S 500 = 0 if there are 0, 2 or more matched sources at 350 pm. 

Third step: deblending and photometry using FastPhot. 

- Use the 250 pm positions obtained in the first step as prior 
positions. 

- Use flux densities at each SPIRE wavelength obtained in the 
second step. 

- Perform simultaneous PSF-fitting and deblending; 

- Assign the measured flux densities that were previously 
missing. 

Table IA.2I presents some statistics following the band¬ 
merging process. 

Appendix B: SPIRE number counts tables 

The measured number counts shown in Fig. 0 are provided in 
tables IB. II IB. 21 and IB. 31 We note that those counts are not cor¬ 
rected for incompleteness or for flux boosting, since here we are 
interested only in the relative distribution of the samples (and the 
quantification of the excess of one sample to another). 


Appendix C: Number counts by type: overdensity 
or lensed sources 

As discussed in Sect. 14.41 Fig. IC.ll renresents the counts for over¬ 
density fields (left) and lensed fields (right). This shows that the 
bright part of the number counts is dominated by single bright 
objects, i.e., the lensed sources. However, these lensed counts 
cannot be used for statistical studies of the lensed sources be¬ 
cause we have only used the surface area of the Planck IN re¬ 
gion. 

Appendix D: Overdensities using AKDE 

We can also use another estimator for the density contrast (>350 
by taking the measured source density from AKDE. If we com¬ 
pute the density minus the median density divided by the median 
density, we obtain the distribution shown in Fig. EH 



AKDE overdensity contrast <S 350 

Figure D.l. Cumulative histogram of the normalized overden¬ 
sity contrast (>350 using the AKDE density estimator. Blue repre¬ 
sents all our SPIRE sources, red represents only redder SPIRE 
sources, defined by S350/S250 > 0.7 and Ssoo/Ssso > 0.6, and 
orange 500 random fields in Lockman. Most of our fields show 
overdensities larger than 5. See Sect. l4.2l for details. 


Appendix E: Gallery of selected sources 

We provide 3-colour images for a few sources drawn from the 
228 of our sample. Each plot shows: 

- a 3-colour image of the SPIRE data, with 250 pm in blue, 
350 pm in green, and 350 pm in red (using the python APLpy 
module: show_rgb); 

- the Planck contour (IN region, i.e., 50% of the peak) as a 
bold white line; 

- the significance of the SPIRE overdensity contrast (starting 
at 2 cr, and then incrementing by 1 cr) as yellow solid lines. 

Figure |E] summarizes the 228 SPIRE sample. 


1 APC, AstroParticule et Cosmologie, Universite Paris Diderot, 
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris 
Cite, 10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 
13, France 


18 
















Planck Collaboration: high-z sources discovered by Planck and confirmed by Herschel 

Table B.l. Differential Euclidean-normalized number counts S 2 ' 5 ^f at 250 j-im. not corrected for incompleteness or flux boosting. 
They are given here only for relative distributions. 


(S 250 ) S 250 range 
[mJy] [mJy] 

IN Planck 
[Jy 1 - 5 so 1 ] 

OUT Planck 
[Jy 1 - 5 sr" 1 ] 

Total Planck 
[Jy 15 sr -1 ] 

Lockman 
[Jy 1 - 5 sr '| 

HLS 

[Jy 1 - 5 sr -1 ] 

11. 

11.- 12. 

— 

— 

— 

— 

— 

12. 

12.- 13. 

- 

- 

- 

- 

20. ± 7. 

13. 

13.- 14. 

- 

- 

- 

- 

103. ± 14. 

16. 

14.- 17. 

- 

- 

- 

- 

271. ± 21. 

19. 

17.- 21. 

235. ± 50. 

476. ± 28. 

445. ± 25. 

- 

447. ± 27. 

25. 

21.- 29. 

9276. ± 343. 

6075. ± 108. 

6493. ± 104. 

4828. ± 122. 

750. ± 38. 

34. 

29.- 40. 

40486. ± 838. 

21789. ±238. 

24227. ± 234. 

29237. ±351. 

3299. ± 94. 

50. 

40.- 59. 

37392. ± 1001. 

18877. ±275. 

21291. ±273. 

23910. ± 395. 

15891. ± 257. 

75. 

59.- 90. 

19583. ± 944. 

10611. ±269. 

11781. ±264. 

13593. ± 388. 

11697. ± 288. 

116. 

90.- 141. 

7545. ± 791. 

4960. ± 248. 

5297. ± 239. 

6494. ± 362. 

6011. ± 278. 

182. 

141.- 223. 

3482. ± 742. 

2587. ± 248. 

2704. ± 236. 

3696. ± 377. 

3127. ± 277. 

290. 

223.- 357. 

935. ± 540. 

1963. ± 303. 

1829. ± 273. 

2730. ± 455. 

1891. ± 303. 

467. 

357.- 576. 

626. ± 626. 

1785. ± 409. 

1634. ± 365. 

1524. ± 482. 

975. ± 308. 

753. 

576.- 931. 

5100. ± 2550. 

573. ± 331. 

1164. ± 440. 

310. ± 310. 

1190. ± 486. 


Table B.2. Differential Euclidean-normalized number counts S 2 5 ^ at 350 yum, not corrected for incompleteness or flux boosting. 
They are given here only for relative distributions. 


(S' 35 o> S 3 5 o range 

IN Planck 

OUT Planck 

Total Planck 

Lockman 

HLS 

[mJy] 

[mJy] 

[Jy 1 - 5 so 1 ] 

[Jy 1 ' 5 so 1 ] 

[Jy 1 - 5 so 1 ] 

[Jy 1 - 5 so 1 ] 

[Jy 1 - 5 so 1 ] 

11 . 

11 .- 12 . 

— 

— 

— 

— 

— 

12 . 

12.- 13. 

- 

- 

- 

- 

- 

13. 

13.- 14. 

- 

- 

- 

- 

5. ± 3. 

16. 

14.- 17. 

- 

- 

- 

- 

39. ± 8 . 

19. 

17.- 21. 

107. ± 34. 

162. ± 16. 

155. ± 15. 

- 

148. ± 16. 

25. 

21.- 29. 

7931. ± 317. 

4114. ± 88 . 

4611. ± 87. 

648. ± 45. 

415. ± 29. 

34. 

29.- 40. 

32868. ± 755. 

13424. ± 187. 

15959. ± 190. 

19667. ± 288. 

3661. ± 99. 

50. 

40.- 59. 

31093. ±913. 

10831. ±209. 

13473. ±217. 

16222. ± 325. 

8986. ± 194. 

75. 

59.- 90. 

16805. ± 875. 

4650. ± 178. 

6235. ± 192. 

6880. ± 276. 

4938. ± 187. 

116. 

90.- 141. 

3731. ± 556. 

1305. ± 127. 

1622. ± 132. 

1936. ± 198. 

1509. ± 140. 

182. 

141.- 223. 
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Figure C.l. Differential Euclidean-normalized number counts, S 25 dN/dS , for various data sets at 350 yum, to illustrate the difference 
between the overdensity fields (left) and the lensed fields (right). This is the same as Fig. [5] except that we have split the fields into 
only overdensity fields and only lensed fields. Lensed sources thus dominate the statistics at large flux densities, while overdensities 
dominate at smaller flux densities. See Sect l4.4l for details. 
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Table B.3. Differential Euclidean-normalized number counts S 2 ' 5 ^f at 500 j-im. not corrected for incompleteness or flux boosting. 
They are given here only for relative distributions. 


<S 5 00 > 5 500 range 

IN Planck 

OUT Planck Total Planck Lockman 

HLS 

[mJy] 

[mJy] 

[Jy 1 - 5 sr -1 ] 

[Jy 1 - 5 sr ‘] 

[Jy 1 - 5 sr -1 ] 

[Jy 1 - 5 sr- 1 ] 

[Jy 1 - 5 sr -1 ] 

11. 

11.- 12. 

_ 

— 

— 

— 

— 

12. 

12.- 13. 

- 

- 

- 

- 

- 

13. 

13.- 14. 

- 

- 

- 

- 

- 

16. 

14.- 17. 

- 

- 

- 

- 

3. ± 2. 

19. 

17.- 21. 

- 

5. ± 3. 

4. ± 2. 

- 

42. ± 8. 

25. 

21.- 29. 

685. ± 93. 

346. ± 26. 

390. ± 25. 

- 

138. ± 17. 

34. 

29.- 40. 

13397. ± 482. 

3474. ± 95. 

4768. ± 104. 

4614. ± 140. 

437. ± 34. 

50. 

40.- 59. 

15600. ± 647. 

2890. ± 108. 

4547. ± 126. 

5692. ± 193. 

1964. ± 90. 

75. 

59.- 90. 

5875. ± 517. 

751. ± 72. 

1419. ± 92. 

1473. ± 128. 

935. ± 81. 

116. 

90.- 141. 

1078. ± 299. 

199. ± 50. 

314. ± 58. 

484. ± 99. 

232. ± 55. 

182. 

141.- 223. 

791. ± 354. 

142. ± 58. 

227. ± 68. 

77. ± 54. 

172. ± 65. 

290. 

223.- 357. 

1559. ± 697. 

140. ± 81. 

325. ± 115. 

76. ± 76. 

97. ± 69. 

467. 

357.- 576. 

2506. ± 1253. 

94. ± 94. 

408. ± 183. 

- 

- 

753. 

576.- 931. 

3825. ± 2208. 

- 

499. ± 288. 

- 

- 
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Figure E.l. Representative Planck targets, showing 3-colour SPIRE images: blue, 250/./m; green, 350/rm; and red, 500/mi. The 
white contours show the Planck IN region, while the yellow contours are the significance of the overdensity of 350 /mi sources, 
plotted 2 cr, 3 cr, 4 <x, etc. 
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FigureE.2. Representative Planck targets. Same as Fig. IE.ll 
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Figure E.3. Representative Planck targets. Same as Fig. IQ 
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Figure E.4. Representative Planck targets. Same caption as Fig. IE.II 

- I’Wt' in, Cncn * ftws - S00*m 

iff 

I 

~W 

Figure E.5. Representative Planck targets. S 



19 Department of Astrophysics/IMAPP, Radboud University 
Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 

20 Department of Physics & Astronomy, University of British 
Columbia, 6224 Agricultural Road, Vancouver, British Columbia, 
Canada 

21 Department of Physics and Astronomy, University College 
London, London WC1E 6BT, U.K. 

22 Department of Physics, Florida State University, Keen Physics 
Building, 77 Chieftan Way, Tallahassee, Florida, U.S.A. 

23 Department of Physics, Gustaf Hallstromin katu 2a, University of 
Helsinki, Helsinki, Finland 

24 Department of Physics, Princeton University, Princeton, New 
Jersey, U.S.A. 

25 Department of Physics, University of California, Santa Barbara, 
California, U.S.A. 

26 Department of Physics, University of Illinois at 
Urbana-Champaign, 1110 West Green Street, Urbana, Illinois, 
U.S.A. 


27 Dipartimento di Fisica e Astronomia G. Galilei, Universita degli 
Studi di Padova, via Marzolo 8, 35131 Padova, Italy 

2X Dipartimento di Fisica e Scienze della Terra, Universita di Ferrara, 
Via Saragat 1, 44122 Ferrara, Italy 

29 Dipartimento di Fisica, Universita La Sapienza, P. le A. Moro 2, 
Roma, Italy 

30 Dipartimento di Fisica, Universita degli Studi di Milano, Via 
Celoria, 16, Milano, Italy 

11 Dipartimento di Fisica, Universita degli Studi di Trieste, via A. 
Valerio 2, Trieste, Italy 

’ 2 Dipartimento di Fisica, Universita di Roma Tor Vergata, Via della 
Ricerca Scientifica, 1, Roma, Italy 

33 Discovery Center, Niels Bohr Institute, Blegdamsvej 17, 
Copenhagen, Denmark 

34 European Southern Observatory, ESO Vitacura, Alonso de Cordova 
3107, Vitacura, Casilla 19001, Santiago, Chile 


22 






Planck Collaboration: high-z sources discovered by Planck and confirmed by Herschel 


hlgp - JSOi-nv (ire*n - JStyrfn, ftirij - SOO|jn 



Figure E.6. Representative Planck targets. Same caption as Fig. IE.II 



Figure E.7. Representative Planck targets. Same caption as Fig. m 
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Figure E.8. Mosaic showing the 228 SPIRE fields. 
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